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ABSTRACT 

The vertical profiles of disc galaxies are built by the material trapped around stable 
periodic orbits, which form their “skeletons”. According to this, the knowledge of the 
stability of the main families of periodic orbits in appropriate 3D models, can predict 
possible morphologies for edge-on disc galaxies. In a pilot survey we compare the 
orbital structures which lead to the appearance of “peanuts” and “X”-like features 
with the edge-on profiles of three disc galaxies (IC 2531, NGC 4013 and UGC 2048). 
The subtraction from the images of a model representing the axisymmetric component 
of the galaxies reveals the contribution of the non-axisymmetric terms. We find a direct 
correspondence between the orbital profiles of 3D bars in models and the observed 
main morphological features of the residuals. We also apply a simple unsharp masking 
technique in order to study the sharpest features of the images. Our basic conclusion 
is that the morphology of the boxy “bulges” of these galaxies can be explained by 
considering disc material trapped around stable 3D periodic orbits. In most models 
these building-blocks periodic orbits are bifurcated from the planar central family of a 
non-axisymmetric component, usually a bar, at low order vertical resonances. In such 
a case the boxy “bulges” are parts of bars seen edge-on. For the three galaxies we 
study the families associated with the “peanut” or “X”-shape morphology are most 
probably bifurcations at the vertical 2/1 or 4/1 resonance. 

Key words: galaxies: kinematics and dynamics; ISM: dust, extinction - Infrared: 
galaxies - Infrared: ISM 


1 INTRODUCTION 

In Patsis & Grosb0l (1996), Patsis, Athanassoula, Grosb0l 
et al. (2002a), and Patsis, Skokos & Athanassoula (2002b, 
hereafter PSA), have been presented boxy orbital profiles 
of 3D, time-independent, models of disc galaxies. They have 
been built by combinations of stable, 3D, orbits belonging in 
most cases to families of periodic orbits bifurcated from the 
planar xl family (see e.g. Contopoulos 2002, or Contopoulos 
& Grosb0l 1989) at the vertical n/1 resonances. The most 
efficient dynamical mechanism for introducing vertical res¬ 
onances in a system is the presence of non-axisymmetric 
components, especially bars. 

Vertical profiles of 3D bars have been constructed by 
PSA, based on the orbital analysis of 3D Ferrers bars (Pfen- 
niger 1984, 1985; Skokos, Patsis V Athanassoula 2002a,b). 
A profile of a single family includes stable orbits in a range 
of “energies” (Jacobi constants - Ej). In order to construct a 
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family profile one needs to know its stable orbits. Out of this 
library of stable periodic orbits any subgroup that helps in 
matching the observed morphological structure can be cho¬ 
sen. Every family has its own characteristic orbital profile, 
and the profiles of the same family in several models are 
similar. Combinations of the profiles of several families of a 
model comprise a model profile. Particular useful represen¬ 
tations of the orbital profiles that allow comparisons with 
real galaxies, or with snapshots of V-body simulations, are 
their weighted images (PSA). The advantage of these images 
is that they show the relative importance of every orbit in 
a profile. They also show the details of the morphology in 
locations where more than one orbit contribute. The images 
are composed by periodic orbits weighted by the mean den¬ 
sity of the model at the points visited by the orbit. Along 
each orbit one picks points at equal time steps. The den¬ 
sity of the model is calculated at each of these points and 
the mean density is taken as the weight of the orbit. Hav¬ 
ing constructed images for every orbit (normalized over its 
total intensity), we can combine them to built a profile for 
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a family of orbits. For the profile of a family are used orbits 
equally spaced in their mean radius. 

The basic conclusions of the orbital analysis relevant to 
the present study are: 

(i) The morphology supported by the composite orbital 
profile of a family may differ from the morphology of the in¬ 
dividual periodic orbits calculated at a particular Ej . Thus, 
in order to get the backbone of a structure to be compared 
with the morphology of a galaxy, one needs to know the 
evolution of the shape of the stable orbits of a family as a 
function of Ej. 

(ii) The radial extent of the profile of an orbital 3D family 
is usually confined within a radius corresponding to a certain 
Ej value. Orbits with a Jacobi constant larger than this 
Ej^ increase their size by increasing practically only in the 
vertical dimension, and this leads to models with stair-type 
edge-on profiles (for more details see Patsis et al. 2002a). 

To these two, one should add that it is the presence of 
the vertical resonances and not the detailed kind of pertur¬ 
bation in the models that shapes the boxiness of the orbital 
profiles. However, “X”-like features as these we discuss here 
are typical of strong bar components. 

Recently, unsharp masking techniques applied to im¬ 
ages of galaxies (Aronica, Athanassoula, Bureau et al. 
2003; Aronica, Bureau, Athanassoula et al. 2004; Bureau, 
Athanassoula, Chung et al. 2004) as well as to images of 
snapshots of A'-body simulations (Athanassoula 2005a,b) 
have shown excellent agreement between the image mor¬ 
phologies and what is predicted by the orbital theory in 
PSA. In particular, Aronica et al. (2003), compared the im¬ 
age of ESO 597-036 after unsharp masking with a model in 
PSA. They found besides a conspicuous “X”-shape feature 
in the central part, surface brightness enhancements along 
the equatorial plane of the galaxy. Both features have their 
counterparts in the orbital models and can be explained 
by material trapped around stable periodic orbits. Simi¬ 
lar features are indicated by Aronica et al. (2004) also for 
ESO-443G042 and by Bureau et al. (2004) for the case of 
NGC 128. The comparison of features between orbital and 
A-body models on the other hand show agreement in even 
finer details. In both papers by Athanassoula (2005a,b) one 
can find in the edge-on views of the models besides “X”-like 
features, density enhancements on the equatorial plane and 
features like “parentheses” (see e.g. Eig. 6 in Athanassoula 
2005a). This is a strong indication that a large percentage of 
the material in the A-body simulation follows orbits around 
the families bifurcated at the vertical n/1 resonances with 
small n. 

In the present paper we apply two image processing 
techniques on the images of three edge-on disc galaxies with 
rather boxy “bulges” in order to detect structures similar 
to those predicted by the orbital theory. We focus our at¬ 
tention to the structures observed in the central regions of 
IG 2531, NGG 4013 and UGG 2048 (NGG 973). Eirst we sub¬ 
tract from the I-band images of the galaxies an axisymmetric 
model that was found by Xilouris, Kylafis, Papamastorakis 
et al. (1997) to describe best the smooth distribution of stars 
and dust in these objects. The models are used to isolate the 
non-axisymmetric term in the profiles of the galaxies. This 
term is expected to reflect in a straightforward way the pres¬ 
ence of vertical resonances (PSA). A second technique we 


apply is just a gaussian filtering (unsharp masking) on the 
images. In Section |2| we give information about the observa¬ 
tions of the galaxies, in Section El we describe the different 
image processing techniques, in Sectional we give the results 
of the image processing we performed and finally we discuss 
our results and present our conclusions in Section El 


2 OBSERVATIONS 

The observations and the data reduction of the galaxies an¬ 
alyzed in this study are presented in Xilouris, Kylafis, Papa¬ 
mastorakis et al. (1997), and Xilouris, Byun, Kylafis et al. 
(1999). We briefly repeat here both. 

Observations of the galaxies UGG 2048 and NGG 4013, 
were made at Skinakas observatory in Grete, using the 1.3 
m telescope, where a Thomson 1024 x 1024 GGD camera 
is installed at the prime focus of the //7.7 Ritchey-Gretien 
telescope. The 19/im pixels of this camera correspond to 0.39 
arcseconds on the sky giving a total field of view 6.7' x 6.7'. 
IG 2531 was observed with the 1 m Australian National Uni¬ 
versity telescope (ANU) at the Siding Spring Observatory. 
The GGD camera in this case is an EEV 576 x 380 giving 
a pixel size of 0.56 arcseconds at the //8 Gassegrain focus. 
The I passband comparable to that of the Gousin’s photo¬ 
metric system is used in both cases. The exposure times are 
20, 25 and 25 minutes for UGG 2048, NGG 4013 and IG 2531 
respectively. 

In Table 1 we summarize the basic observational proper¬ 
ties of the galaxies as well as some characteristic parameters 
describing the axisymmetric stellar and dust components 
(see Section El . 

Out of the samples of the galaxies in Xilouris et al. 
(1997) and Xilouris et al. (1999), for which models for the 
axisymmetric light and dust distribution exist, we have se¬ 
lected three of them with boxy profiles to apply our image 
processing techniques. 


3 DETERMINING THE UNDERLYING 

STRUCTURE OF THE GALAXY 

Disc galaxies are not axisymmetric systems. Their surface 
brightness can be described as the sum of an axisymmetric 
and a perturbing term. The latter refers to the presence of 
spirals or bars. The axisymmetric component is the most 
important and can be described by simple smooth functions 
(like axisymmetric exponential discs for the stars and the 
dust in the plane of the galaxy and a de Vaucouleurs 
law for the bulge). Besides the global perturbations we have 
clumps and star forming regions that cause deviations from 
smoothness locally. In the present paper we are interested 
in the global non-axisymmetric components. 

Studying the smooth distributions is not only instruc¬ 
tive (one gets a very good description of the galaxy with as 
few parameters as possible) but can also be a very helpful 
tool to uncover features due to non-axisymmetric compo¬ 
nents (e.g. bars). This is done by subtracting the axisymmet¬ 
ric model from the image of the galaxy. If the axisymmetric 
model is able to describe well the large-scale structure, then, 
with this method, what it left, i.e. the residuals, will indicate 
features that mainly have to do with either the non-smooth 
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Table 1. Observational data for the galaxies together with a set of parameters describing the axisymmetric stellar and dust components 
(see the text in Section l^for a detailed description of these parameters). 


Name 

Pixel size 

// 

Exp. time 
(min) 

D 

(Mpc) 

Inclination 

(degrees) 

he 

(kpc) 

hd 

(kpc) 

Zs 

(kpc) 

Zd 

(kpc) 

Re 

(kpc) 

Is 

(mag/o") 

h 

(mag/n”) 

UGC 2048 

0.39 

20 

63 

89.6 

11.0 

16.5 

1.0 

0.57 

2.4 

18.4 

9.1 

NGC 4013 

0.39 

25 

12 

89.6 

1.8 

2.6 

0.2 

0.10 

1.6 

17.3 

10.5 

IG 2531 

0.56 

25 

22 

89.7 

5.0 

8.4 

0.4 

0.20 

1.6 

18.2 

11.0 


distribution of the dust (seen as absorption features) or the 
non-axisymmetric component of the stars (seen as excess of 
starlight). Assuming a constant mass-to-light ratio (M/L), 
this light excess will indicate regions of enhanced stellar den¬ 
sity and will be comparable with the orbital models. Good 
agreement will indicate that large percentages of stars are 
trapped around stable periodic orbits. 

The model that we use is described in great detail in 
Xilouris et al. (1997). The stellar emissivity (luminosity per 
unit volume) that we use consists of an exponential (in both 
radial and vertical directions) disc and a bulge described by 
the law, namely 

\ hs Zs J 

+ u exp(-7.67B^/^)B“^/^ (1) 


with hs and Zs being the scalelength and scaleheight of the 
disc and 


B = 


Re 


( 2 ) 


with Re being the effective radius of the bulge and a and 
b the semi-major and semi-minor axis respectively. Here Lg 
and Lb are the normalization constants for the stellar emis¬ 
sivity of the disc (Lg) and the bulge (Lb). The central value 
for the surface brightness of the disc and the bulge, if the 
model galaxy is seen edge-on and there is no dust, are given 
by Ig = 2Lghg and R = 5.12LbRe- 

For the extinction coefficient we use a double exponen¬ 
tial law, namely 


k\{R,z) = exp 




( 3 ) 


where Hi\ is the extinction coefficient at wavelength A at the 
center of the disc and hd and Zd are the scalelength and 
scaleheight respectively of the dust. 

The radiative transfer model that we have used is the 
one described by Kylahs & Bahcah (1987) (see also Xilouris 
et al. 1997). The results are summarized in Table 1, where 
we give the values of the most important parameters in the 
best httings. 

The unsharp masking has been done by means of stan¬ 
dard ESO-MIDAS commands. First we apply on the image 
a gaussian hlter with radius of 5 pixels around the central 
pixel and with mean and sigma values of 5 and 2 pixels re¬ 
spectively in both directions. With this procedure a blurred 
image of the galaxy is created. Then, after subtracting the 
hltered from the original image the residual indicates the 
presence of the sharpest features on the images. 


4 RESULTS 

In Fig. 1 we show the residual maps created by subtracting 
the model images of Xilouris et al. (1997, 1999) from the ob¬ 
served images of the three galaxies (NGC 4013 (a), IG 2531 
(b) and UGG 2048 (c)). The images have been rotated so 
that the major axis of the galaxies is along the large dimen¬ 
sion of the frame. The difference between the models and 
the observations are relatively small, indicating the “pertur¬ 
bation” character of the non-axisymmetric components (on 
average 15% - Xilouris et al. 1997, 1999). However, this pro¬ 
cedure offers a straightforward tracer exactly of this excess 
light we are interested in. The main characteristic feature 
is an “X”-shaped structure revealed in the central parts of 
the galaxies. By comparing the residual with the full images 
of the galaxies we observe that the “X” features are in the 
region, where we have the boxy bulges. 

The images after unsharp masking show sharp features 
that resemble the wings of an “X”. This can be seen in Fig. 2 
for the same galaxies. Gharacteristic contours underline the 
main structure observed. In this case we can better distin¬ 
guish the highest contrast features. Again they are located 
at the region occupied by the boxy bulges. They should be 
compared to features expected in the central parts of 3D 
bars viewed almost, edge-on. The “X-shape” is encountered 
only in models for bars, while boxy bulges in general can be 
found also in non-barred models (Patsis et al. 2002a). 


5 DISCUSSION AND CONCLUSIONS 

Both image processing techniques revealed a kind of “X” 
structure in the central components of the three edge-on 
galaxies with boxy prohles we studied. In Fig. 1, where we 
give the residuals after subtracting a Xilouris et al. (1997) 
model we see the overall shape of the non-axisymmetric com¬ 
ponent. The unsharp masking in Fig. 2 shows more clearly 
the high intensity features of the “bulges”. They emerge out 
of the equatorial plane as distinct branches. 

In most cases of the models in PSA an “X” is sup¬ 
ported by material trapped around xlvl orbits (Skokos et 
al. 2002a), i.e. by stable 3D orbits introduced in the system 
at the vertical 2/1 resonance. However, there are other 3D 
families as well, which could support straight line segments 
emerging out of the equatorial plane in the edge-on projec¬ 
tions of the orbital models. Such is the case of xlv5, accord¬ 
ing to the nomenclature of Skokos et al. (2002a), which is a 
family born at the vertical 4/1 resonance. This family has 
been associated with boxy central components already by 
Pfenniger (1985). A typical orbital prohle is given in Fig. 7a 
in PSA. 

In models (orbital or A-body), a visual difference can 
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Figure 1. Residual images for NGC 4013 (a), IC 2531 (b) and UGC 2048 (c). From the original images it is subtracted a best fitting 
axisymmetric profile according to the models of Xilouris et al. (1997). An “X”-shaped feature embedded in the boxy “bulge” is conspicuous 
in all three cases. Golor coding can be seen below (a). 



Figure 2. The result of unsharp masking on the images of NGG 4013 (a), IG 2531 (b) and UGG 2048 (c). “X” shape features similar 
to the structures foreseen by the orbital models can be found in all three cases. Arrows indicate breaks of the “X” wings. Golor coding 
can be seen below (a). 


be found if we compare the side-on views. The distance be¬ 
tween the wings of the “X” feature in the xlv5 case is larger 
than the distance of the corresponding wings in xlvl pro¬ 
files. The difference can be seen in Fig. 3, where we give in 
(a) a side-on typical xlvl profile and in (b) a profile domi¬ 


nated by the presence of the xlv5 family. They correspond 
to models D and B in PSA. As we mentioned in the intro¬ 
duction, there is a certain value of Ej^ beyond which the 
orbits of a family grow practically only in the z direction. 
So, for each family, there is a maximum radius on the equa- 
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Figure 3. (a) A typical xlvl side-on profile based on orbits of 
model D in PSA. Arrows indicate the inward bending of the “X” 
wings at large distances above the equatorial plane, (b) A typical 
side-on xlv5 profile corresponding to model B in PSA. It extends 
to larger radii than the profile in (a). 

torial plane, within which the projections of its orbits are 
confined. In the side-on views we can measure the length of 
the projection along the major axis of the bar. In the two 
specific examples we give here, the projection of the xlvl 
orbits on the major axis of the bar in Fig. 3a corresponds 
roughly to a length 45% of the longest bar supporting or¬ 
bits, i.e. the length of the bar in our orbital model. On the 
other hand in Fig. 3b, this percentage reaches almost 90%. 
Unfortunately this criterion cannot in general apply on the 
images of real galaxies since we are missing the essential in¬ 
formation about the orientation of the boxy structure with 
respect to the line-of-sight. Thus, at the level of the current 
analysis, we cannot point to a single vertical resonance and 
attribute to it the observed morphological feature. A galaxy 
may lack a vertical 2/1 (or even 3/1) resonance, in which 
case material may be trapped by stable orbits bifurcated 
at the vertical 4/1 resonance. The corresponding profile in 
such a case will look like what we see in Fig. 3b. It con¬ 
sists of stable orbits of the following families: (a) xlv5, (b) 
a bifurcation of it, and (c) z3.1s (PSA). 

Another point that has to be mentioned, is the inwards 
bending or break of the wings of “X” close to their maximum 
height above the equatorial plane. According to the orbital 
models this can be explained by the fact that the orbits 
of a family increase their size practically only in the vertical 
direction beyond a certain Ej value. In Fig. 3a we point with 
arrows to these breaks at the xlvl profile. However, this is 
expected to happen to the profiles of other families as well 
if they are populated with orbits with high energies (Ej) 
(see a characteristic case in PSA, Fig. 18). This tendency 
can be observed at the wings of the “X” features in Fig. 2a 
(NGC 4013) and Fig. 2c (UGC 2048). We indicate these 
breaks with arrows. 


The main conclusion of the present study is that the 
careful subtraction of axisymmetric components from the 
profiles of the three edge-on disk galaxies (NGG 4013, 
IG 2531 and UGG 2048) reveals an “X”-like morphology, 
which indicates the presence of a bar. With this method we 
could isolate the light coming from the non-axisymmetric 
components. The similarity of these components with the 
morphology of the orbital prohles is a strong indication that 
the boxy “bulges” of these three galaxies are parts of bar 
structures observed edge-on. If the “X”-hke structure of the 
residuals is due to stars trapped around stable periodic or¬ 
bits, then we expect the presence of sharp features (the 
wings of the “X”) inside the bar region. By applying the 
unsharp masking on the images we hnd, in the foreseen by 
the orbital models regions, sharp features. The two methods 
are complementary. The unsharp masking does not give any 
information about the overall shape of the non-axisymmetric 
component. It shows however that in the regions we see 
the residuals we have the expected sharp features. These 
features show a nice correspondence with the dense parts 
of published orbital models (see PSA, Fig.la, 3b, 7a, 9a). 
These techniques, and especially the study of the residuals 
after subtracting axisymmetric models, should be applied in 
a large sample of galaxies in order to estimate the frequency 
of these features in the prohles of boxy edge-on galaxies. 
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